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By imposing the constraints of charge symmetry we show that the strangeness magnetic moment of the nucleon 
can be expressed in terms of empirical magnetic moments and ratios of valence quark magnetic moments. The 
latter are determined using modern chiral extrapolation techniques and recent low mass lattice QCD simulations 
of the individual quark contributions to the magnetic moments of the nucleon octet. The result is a precise 
determination of G^/, namely —0.043 ±0.026 hn, which is consistent with the latest experimental measurements. 



1. INTRODUCTION 

There is currently enormous interest in the de- 
termination of the strangeness content of the nu- 
cleon. It is clearly crucial to our understanding 
of QCD to determine precisely the role played by 
heavier, non-valence flavors. On the experimen- 
tal side there have been tremendous advances in 
the ability to measure parity violation in electron 
scattering at the level of 10^^ and new results on 
strangeness in the nucleon have been reported re- 
cently from JLab (HAPPEX) 1 and MIT-Bates 
(SAMPLE) 2 . In the near future we can expect 
even more precise results from the A4 experiment 
at Mainz as well as GO and HAPPEX2 at JLab. 

In contrast, the theoretical situation is some- 
what confused with the predictions of various 
quark models covering an enormous range. Di- 
rect calculations within lattice QCD have not yet 
clarified the situation, with values for G%j ranging 
from -0.28 ±0.10 3 to ±0.05 ±0.06 4 • We take 
a different approach, building on the improve- 
ments in both lattice actions and computer speed 
which have enabled quenched QCD (QQCD) sim- 
ulations of magnetic moments at pion masses as 
low as 0.3-0.4 GeV and on the developments of 
modern chiral extrapolation techniques which al- 
low one to rigorously ensure the model indepen- 
dent constraints of chiral symmetry. Using these 



techniques we determine the ratios of the w-quark 
contribution to the magnetic moment of the phys- 
ical proton to that in the E+ and of the u quark in 
the physical neutron to that in the 'EP. From these 
ratios, experimental data on the octet moments 
and charge symmetry, which is typically satisfied 
at the level of 1% or better j^, we deduce a new 
theoretical value for G^^ which is extremely pre- 
cise - setting a tremendous challenge for the next 
generation of parity violation experiments. 



2. CHARGE SYMMETRY 

An examination of the symmetries manifest 
in the QCD path integral for current matrix el- 
ements reveals various relationships among the 
quark contributions [Hj. The magnetic moment 
of the proton, is extracted from the three-point 
function, where an operator exciting the proton 
from the QCD vacuum is followed by the electro- 
magnetic current, which in turn is followed by an 
operator annihilating the proton back to the QCD 
vacuum. In calculating this three point func- 
tion, one encounters two topologically distinct 
ways of performing the electromagnetic current 
insertion. Figure Q displays skeleton diagrams for 
these two possible insertions (with Euclidean time 
increasing to the right). In full QCD these dia- 
grams incorporate an arbitrary number of gluons 



Figure 1. Diagrams illustrating the two topologi- 
cally different insertions of the current within the 
framework of lattice QCD. These skeleton dia- 
grams for the connected (left) and disconnected 
(right) current insertions may be dressed by an 
arbitrary number of gluons and quark loops. 



and quark loops. The left-hand diagram illus- 
trates the connected insertion of the current to 
one of the "valence" quarks of the baryon. In the 
right-hand diagram the external field produces a 
qq pair which in turn interacts with the valence 
quarks of the baryon via gluons. It is important 
to realize that within the lattice QCD calculation 
of the loop diagram on the right side of Fig. Q] 
there is no anti-symmetrization (Pauli blocking) 
of the quark in the loop with the valence quarks. 
For this reason, in general only the sum of the 
two processes in Fig. Q] is physical. 

Under the assumption of charge symmetry, 
the three-point correlation functions for octet 
baryons leads to the following equalities for elec- 
tromagnetic current matrix elements [H]: 
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= CuuP + cd dP + On 


n 


= eduP + Cu dP + On 


S+ 


= e„ u^ + ess^ + Oe 


E" 


= edu^ + esS^ + 0^ 


^0 


= CsS" + e„ u" + Ob. 


,'— ', 


= Bs s" -1- Cd u" + Ob. 



(1) 

Here, all quantities refer to magnetic moments, so 
for example p and S~ are the physical magnetic 
moments of the proton and S~, respectively. The 
valence w-quark magnetic moment in the proton, 
corresponding to the left-hand side of Fig. Q] is 
denoted vP . We emphasize that charge symmetry 
has been used to replace the d-quark contribution 
in the neutron by u^, d in the S^ by u in the E+ ( 



M^), and so on. Finally, O denotes the total con- 
tribution from quark- loops - i.e., the contribution 
shown on the right-hand side of Fig. ^ The la- 
bels on quark magnetic moments allow for the en- 
vironment sensitivity implicit in the three-point 
function 6 . For example, the three-point func- 
tion for the S+ is the same as that for the proton, 
except that the d is replaced by the somewhat 
heavier s-quark. Hence, the w-quark propagators 
in the S^ are multiplied by an s-quark propaga- 
tor, whereas in the proton they are multiplied by 
a d-quark propagator. The different mass of the 
neighboring quark gives rise to an environment 
sensitivity in the u-quark contributions to observ- 
ables, which means that the naive expectations of 
the constituent quark model vF jv^ = u^ jvr = 1 
may not be satisfied 6 7 8l9ll()lllll'^in] . This 
observation should be contrasted with the usual 
quark model assumption that the quark magnetic 
moment is an intrinsic property, independent of 
the quark's environment. 

In three-flavor QCD, On contains sea-quark- 
loop contributions from u, d and s quarks. In the 
SU(3)-flavor limit {rriu = nid = rris) the charges 
add to zero and hence the sum vanishes. However, 
the heavier strange quark mass leads to a non- 
zero result. By definition 
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and the leading superscript, £, reminds one that 
the contributions are loop contributions. Note 
that, in deriving Eq.©, we have set ^GJ^j — ^Gff , 
corresponding to ?Ti„ — irid [HI- Earlier esti- 
mates of ^i?^ were based on the constituent quark 
model, however a more reliable approach is to 
estimate the loops using a finite range regulator 
|14I15I16) . Since the chiral coefficients for the d 
and s loops in the RHS of Fig. 1 are identical, 
the only difference comes from the mass of the K 
compared with that of the tt. 



With no more than a httle accounting, the 
strange-quark loop contributions to the nucleon 
magnetic moment, G^j may be isolated from ^ 
and (PI in the following two phenomenologically 
useful forms, 
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As we have explained, under the assumption that 
quark magnetic moments are not environment de- 
pendent, these ratios (i.e. u^/u^ and u^ ju") are 
taken to be unity in many quark models. In- 
corporating the experimentally measured baryon 
moments leads to: 
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where all moments are expressed in nuclear mag- 
netons {(J-n)- (Note that the measured magnetic 
moments are all known sufficiently accurately 1171 
that the experimental errors play no role in our 
subsequent analysis.) We stress that these ex- 
pressions for G\i are exact consequences of QCD, 
under the assumption of charge symmetry. Equa- 
tion © provides a particularly favorable case for 
the determination of G^j with minimal depen- 
dence on the valence-quark ratio. 

If one considers the quark model suggestions 
of w"/u" — 1 and ^R^ — 0.65 in ((SJ, one finds 
G\j = —0.81 /iAT, a significant departure from 
the experimental preference of positive values. 

Equating Q and JHJ provides a linear relation- 
ship between vP jiP and u^'' ju^ which must be 
satisfied within QCD. Figure [3 displays this re- 
lationship by the dashed and solid line, the lat- 
ter corresponding to values for which G|,j(0) > 
when ^i?^ is in the anticipated range < 
^i?^ < 1. Since the line does not pass through 
the point (1.0,1.0) corresponding to the simple 
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Figure 2. The consistency relation between 
u^'/u and u"/m" which must be satisfied within 
QCD. The part of the straight line which is 
dashed corresponds to G^^(O) < 0, while the solid 
part of the line has G^^(O) > 0. The standard 
quark model assumption of intrinsic quark mo- 
ments independent of their environment is indi- 
cated by the filled square. The lattice QCD pre- 
diction (after an appropriate chiral extrapolation, 
discussed in following sections) is illustrated by 
the filled circle. 



quark model assumption of universality, the ex- 
perimentally measured baryon moments are sig- 
naling that there must be an environment effect 
exceeding 12% in both ratios or approaching 20% 
or more in at least one of the ratios. More- 
over, a positive value for G^j(O) requires an envi- 
ronment sensitivity exceeding 70% in the u"/m" 
ratio. Hence the experimental suggestion that 
G^(0) > challenges the intrinsic magnetic mo- 
ment concept which is fundamental to the con- 
stituent quark model. 

3. NUMERICAL SIMULATIONS WITH 
FLIC FERMIONS 

Access to leading-edge supercomputing re- 
sources coupled with advances in the formu- 
lation of computationally-inexpensive chirally- 
improved lattice fermion actions enable the nu- 
merical calculation of hadron structure near the 
chiral regime. The numerical simulations of the 



electromagnetic form factors presented here are 
carried out using the Fat Link Irrelevant Clover 
(FLIC) fermion action |18I19| in which the irrele- 
vant operators introduced to remove fermion dou- 
blers and lattice spacing artifacts are constructed 
with smoothed links. These links are created via 
APE smearing [2111; a process that averages a link 
with its nearest transverse neighbors in a gauge 
invariant manner. Iteration of the averaging pro- 
cess generates a "fat" link. 

The use of links in which short-distance fluc- 
tuations have been removed simplifies the deter- 
mination of the coefficients of the improvement 
terms in both the action and its associated con- 
served vector current. Perturbative renormaliza- 
tions are small for smeared links and the mean- 
field improved coefficients used here are sufficient 
to remove 0{a) errors, in the lattice spacing a, 
from the lattice fermion action. The key is that 
both the energy dimension-five Wilson and Clover 
terms [2] are constructed with smooth links, 
while the relevant operators, surviving in the con- 
tinuum limit, are constructed with the original 
untouched links generated via standard Monte 
Carlo techniques. 

FLIC fermions provide a new form of nonper- 
turbative 0{a) improvement |19I22| where near- 
continuum results are obtained at finite lattice 
spacing. Access to the light quark mass regime 
is enabled by the improved chiral properties of 
the lattice fermion action. The magnitude of ad- 
ditive mass renormalizations is suppressed |^ 
which otherwise can lead to singular behavior in 
the propagators as the quarks become light. 

The C'(a)-improved conserved vector current 
|23| is used. Nonperturbative improvement is 
achieved via the FLIC procedure where the terms 
of the Noether current having their origin in the 
irrelevant operators of the fermion action are con- 
structed with mean-field improved APE smeared 
links. The preliminary results presented here are 
from a sample of 255 20^ x 40 mean- field improved 
Luscher-Weisz PH gauge field configurations hav- 
ing a lattice spacing of 0.128 fm as determined by 
the Sommcr scale vq = 0.50 fm. 



4. CHIRAL EXTRAPOLATION 

One of the major challenges at present in con- 
necting lattice calculations of hadronic properties 
with the physical world is that computational lim- 
itations restrict the accessible quark masses to 
values much larger than the physical values. At 
quark masses typical of today's lattice QCD simu- 
lations, one is outside the region where traditional 
dimensionally regulated (DR) chiral perturbation 
theory (xPT) is applicable. Yet one knows that 
for current quark masses near zero there is impor- 
tant non-analytic structure (as a function of the 
quark mass) which must be treated correctly if we 
are to compare with physical hadron properties. 

Our present analysis of the strangeness mag- 
netic form factor has been made possible by a 
significant breakthrough in the regularization of 
the chiral loop contributions to hadron observ- 
ables |14I15| . Through the process of regulating 
loop integrals via a finite-range regulator ^| , the 
chiral expansion is re-summed producing an ex- 
pansion with vastly improved convergence prop- 
erties. 

The key feature of finite-range regularization 
(FRR) is that FRR schemes have an additional 
adjustable regulator parameter which provides an 
opportunity to suppress short distance physics 
from the loop integrals of effective field theory. 
This short-distance physics is otherwise treated 
incorrectly as the naive effective fields do not 
share the properties of QCD at short distances. 

An extensive study of the quark mass depen- 
dence of the nucleon mass in finite-range reg- 
ularized (FRR) xPT UD explored six differ- 
ent regularization and associated renormaliza- 
tion schemes. The smooth regulators (dipole, 
monopole, or Gaussian regulators) provide ex- 
pansions that agree over the range < m^ <~ 
0.8 GeV^ at a level sufficient to predict the 
nucleon mass within 1%. Our focus here is 
to extrapolate FLIC fermion calculations of va- 
lence quark contributions to baryon moments 
(m^, u", u^, u") to the physical mass regime. 
We select the dipole-vertex FRR with A — 0.8 
GeV as determined in Ref. |2S1- This scale was 
found to give the best simultaneous description of 
both quenched and dynamical simulation results. 
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Figure 3. Diagrams providing the leading con- 
tributions to the chiral expansion of baryon mag- 
netic moments (upper diagrams) and their asso- 
ciated quark flows (lower diagrams) in QQCD. In 
the lower diagrams, the quark-photon coupling is 
summed over adjacent quark lines. 



Separation of the valence and sea-quark-loop 
contributions to the meson cloud of full QCD 
hadrons is a non-trivial task. We use the dia- 
grammatic method for evaluating the quenched 
chiral coefficients of leading nonanalytic terms in 
heavy-baryon quenched xPT |26I27| . These re- 
sults are generalized to the FRR approach used 
here. The diagrammatic technique provides a 
transparent means to accomplish the separation 
of magnetic moment contributions in fuU-QCD 
into "direct sea-quark loop" and "valence" con- 
tributions as in the discussion surrounding Fig.^ 
The valence contributions (key to this analysis) 
are obtained by removing the direct-current cou- 
pling to sea-quark loops from the total contribu- 
tions. Upon further removing "indirect sea-quark 
loop" contributions, where a valence quark forms 
a meson composed with a sea-quark loop, one ob- 
tains the "quenched valence" contributions, the 
conventional view of the quenched approxima- 
tion. Isolating a particular quark flavour contri- 
bution only requires setting the electric charge of 
all other quark flavours to zero. 

Figure displays the diagrams providing the 
leading contributions to the chiral expansion of 
baryon magnetic moments (upper diagrams) and 
their associated quark flows in quenched QCD 



(QQCD). The associated chiral expansion for the 
proton magnetic moment, /ip, has the form 

Up = Oq -|-/Xj,Xj)'-f>)'("^7r, A) +X7r/7r(m^,A) 

+XKlK{mK,K)+4ml + a^ml. (9) 
where / denotes a loop integral defined by 
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(11) 



/^/(m^,A) = -/ dk -u^{k,A) (12) 

Jo (fc2+m2)2 

where A^^v is the physical mass-splitting be- 
tween the N and S or A. The coefficients, x, 
denote the model-independent coefficients of the 
LNA term for tt and K mesons !^. We take 



m|f = m'-j^' " + -171^ 



(13) 



where the physical values may be used to define 
The m^ term of Eq. J^ allows for some 



4^ 



curvature associated with the Dirac moment of 
the baryon ex i/m'^ for moderately large quark 
masses. 

The un-renormalized coefficients of the ana- 
lytic terms of the FRR expansion are regulator- 
parameter dependent. This is emphasized by the 
superscript A on the coefficients oq, 0,2 and 04. 
The large mjr behavior of the loop integrals of 
Eq. (|5J and the residual expansion are remark- 
ably different. Whereas the residual expansion 
will encounter a power divergence, the FRR loop 
integrals will tend to zero as some power of A/m^, 
as ?7T,7r becomes large. Thus, the A dependence of 
flQ, a2 and 04 provides an opportunity to govern 
the convergence properties of the residual expan- 
sion and thus the FRR chiral expansion. 

Lattice QCD has now provided extensive 
model-independent information on the moder- 
ate to large ttItt dependence of hadron observ- 
ables. In particular, hadron masses are observed 



to be smooth almost linear functions of 



for 



quark masses similar to the strange quark mass. 
The magnetic moments presented here are also 
smooth, taking on a Dirac moment dependence 



as the quark mass becomes large. This indicates 
that it should be possible to tune the regulator- 
range parameter, A, such that the coefficients a4, 
and higher are truly small. In this case the con- 
vergence properties of the residual expansion, and 
the loop expansion are excellent and their trun- 
cation benign. Indeed this hypothesis was con- 
firmed for the nucleon mass in Ref. |14j . 

It is important to note that the optimal A is one 
that optimizes the convergence properties of the 
residual expansion. A poor choice for A will move 
strength in terms proportional to 1 , rn^ , m* , . . . 
from the loop integrals into the residual expan- 
sion, changing the convergence properties of the 
FRR expansion. We emphasize that the optimal 
A is not selected to approximate the higher-order 
terms of the chiral expansion. If an infinite num- 
ber of parameters were being approximated by 
one, one would have an uncontrolled error. For- 
tunately such ambitions are not necessary. The 
lattice QCD results indicate that these higher or- 
der terms of the chiral expansion largely sum to 
zero. Moreover, the findings of Ref. J^ indicate 
that the details of exactly how each of the terms 
enters the vanishing sum are largely irrelevant. 

Perhaps it is also worth emphasizing that the 
choice of A can have no effect on the convergence 
properties of the renormalized chiral expansion. 
There, contributions from the residual expansion 
and the FRR loop integrals are combined, remov- 
ing any A dependence from the chiral expansion. 
This is the process of renormalization. 

Donoghue et al. [2H| describe the convergence 
problems of DR as due to incorrect short-distance 
physics in loop integrals which must be removed 
by large and opposing analytic-term contribu- 
tions. The introduction of an ultraviolet cut off in 
FRR schemes prevents the error from being made 
in the first place and provides the opportunity to 
have excellent convergence properties in the resid- 
ual expansion. In other words, the higher-order 
analytic terms of the FRR residual expansion do 
not need to be large, as incorrect short distance 
physics has been suppressed from the loop inte- 
grals. The beauty of the FRR procedure is that 
upon renormalization, very large coefficients can 
be recovered [T^ . 

Figure H illustrates a fit of FRR QxPT to the 
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Figure 4. The contribution of a single u quark 
normalized to unit charge to the magnetic mo- 
ment of the proton. Lattice simulation results 
(square symbols for mj > 0.05 GeV) are extrap- 
olated to the physical point (vertical dashed line) 
in finite-volume QQCD with a discrete momen- 
tum sum (solid curve, square symbol), and in 
infinite-volume QQCD (long-dashed curve, dia- 
mond symbol) . Estimates of the valence u quark 
contribution in full QCD (dot-dash curve, trian- 
gular symbol) and the full u quark sector contri- 
bution in full QCD (fine-dash curve, circular sym- 
bol) normalized to a single quark unit charge are 
also illustrated. Extrapolated values indicated 
by symbols at the physical pion mass (vertical 
dashed line) are offset for clarity. 



FLIC fermion lattice results (solid curve) , in a cal- 
culation where only the discrete momenta allowed 
in the finite volume of the lattice are summed in 
performing the loop integral. The long-dashed 
curve that also runs through the lattice results 
illustrates the case when the discrete momentum 
sum is replaced by the infinite- volume, continu- 
ous momentum integral. For all but the lightest 
quark mass, finite volume effects are negligible. 

The coefficients of the residual expansion, 
ag , 02 , a^ , show excellent signs of convergence 
with values 1.52(14), -1.23(48), and 0.69(42) in 
appropriate powers of GeV, respectively. 

Incorporating baryon mass splittings into the 
kaon loop contributions is essential. For example 
the contribution of S ^ NK is nearly doubled 




Figure 5. Correcting QxPT (upper diagrams) to 
full QCD (lower diagrams). Quenched negative- 
metric rj' contributions must be removed, while 
the chiral coefficients of tt and K loops must be 
adjusted to account for the new diagram where 
a valence quark couples to the photon in a me- 
son constructed from a sea-quark loop. Photon 
couplings to the anti-quark in the bottom-right 
diagram are also included in the valence quark 
contributions to the nucleon moment of full QCD. 
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Figure 6. The contribution of a single u quark 
normalized to unit charge to the magnetic mo- 
ment of E+ . Curves and symbols are as described 
in Fig. El 



when the S — iV mass splitting is taken into ac- 
count. 

Figure diagrammatically illustrates the con- 
siderations in correcting the quenched w-quark 
contribution to the valence u-quark contribution 
in full QCD. The removal of quenched rj' contri- 
butions and the appropriate adjustment of tt and 
K loop coefficients as detailed in Ref. 26 27^ pro- 
vides the dot-dash curve of Fig. 0| This is our 
best estimate of the valence u-quark contribution 
(the connected insertion of the current) to the 
proton magnetic moment of full QCD. The graph 
corresponds to a single u-quark contribution nor- 
malized to unit charge. It will be interesting to 
confront this curve with full QCD simulation re- 
sults in the future. 

Finally we also include the disconnected inser- 
tion of the current in order to estimate the total 
contribution of the u-quark sector to the proton 
magnetic moment. This is represented by the fine 
dashed curve in Fig. ^ where the sector contribu- 
tion is normalized to a single quark of unit charge. 
Figures El U\ and |H1 display similar results for the 
E+, n and E°. 
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Figure 7. The contribution of the u quark nor- 
malized to unit charge to the magnetic moment of 
the neutron. Curves and symbols are as described 
in Fig. 21 
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Figure 8. The contribution of the u quark nor- 
mahzed to unit charge to the magnetic moment 
of "EP . Curves and symbols are as described in 
Fig.H 
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Figure 9. Comparison of the FRR^PT extrapo- 
lated lattice QCD simulation results (square sym- 
bols) with experimentally measured baryon mag- 
netic moments (circular symbols). 



From these chiral extrapolations, we estimate 
the ratio of the valence (connected) u-quark con- 
tribution to the proton to that in S+. Simi- 
larly, the ratio of the u-quark contribution to the 
neutron relative to that in S° is calculated. Er- 
ror bars are estimated from a third-order single- 
elimination jackknife analysis of the correlated ra- 
tios. Our final result of 



1.07 ±0.04 and 



1.35 ±0.15 (14) 



is plotted in Fig. |21 This result leaves little doubt 
that G%j is negative. The fact that this point lies 
exactly on the constraint curve is highly nontriv- 
ial, and provides a robust check on the validity of 
the analysis techniques presented here. 

As a further check we compare the lattice QCD 
predictions of the baryon magnetic moments con- 
structed from extrapolations of the individual 
quark sectors in Fig. El The results display an un- 
precedented agreement with experiment. There 
is a hint of a small systematic error that sup- 
presses the singly represented quark contribution 
to both S baryons and the neutron. We have 
identified contributions from decuplet baryons as 
key to resolving this minor discrepancy. However, 
the ratio of m"/u" does not appear to be affected 
by this contribution, as evidenced by the excel- 
lent agreement between lattice and experiment 



for the u quark in S and the degree of consis- 
tency between the lattice extrapolations and the 
constraint curve of Fig. |21 The fact that system- 
atic errors will cancel in taking the ratios of spe- 
cific quark contributions was the main motivation 
for introducing the ratios in the first place. 

While G%j is most certainly negative, it re- 
mains to set the magnitude. This requires an 
estimate of the strange to light sea-quark loop 
contributions, ^R^- Earlier estimates oi^R^ were 
based on the constituent quark model, however 
a more reliable approach is to estimate the loops 
using a finite range regulator |14ll5j . Using a 
dipole-vertex regulator, a direct calculation pro- 
vides 

^i?^ = ^ =0.160 ±0.060, (15) 

where the range includes the variation of the 
dipole mass parameter between 0.6 and 1.0 GeV 
as well as the possible variation of the mass of 
the strange baryon accompanying the K^ loop. 
Using this value in Eq. (jBJ with the results of 
Eq. H14|) provides the final precise determination 
of 



GIj = -0.043 ± 0.026 UN 



(16) 

for the strange quark contribution to the mag- 
netic moment of the nucleon. This extremely pre- 



cise value sets a tremendous challenge for the next 
generation of parity violation experiments. 
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